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Abstract 



m 

> 

| The paper deals with the integral equation approach to steady kinematic 

dynamo models in finite domains based on Biot-Savart's law. The role of the 

o ' 

electric potential at the boundary is worked out explicitly. As an example, 
Q-f a modified version of the simple spherical a-effect dynamo model proposed 

6 

by Krause and Steenbeck is considered in which the a-coefficient is no longer 
constant but may vary with the radial coordinate. In particular, the results 
for the original model are re-derived. Possible applications of this integral 
equation approach for numerical simulations of dynamos in arbitrary geome- 
try and for an "inverse dynamo theory" are sketched. 



I. INTRODUCTION 

For decades, theoretical and numerical work on dynamos has been done for the main part 
in terms of differential equation systems. Countless computer simulations, in particular for 
spherical geometry, have led to a good understanding of dynamos, at least on the kinematic 
level. 
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In a few papers the steady state of dynamos in an infinitely extended fluid with constant 
conductivity has been investigated on the basis of integral equations; see Gailitis (1967), 
Gailitis (1970), Gailitis and Freiberg (1974), Dobler and Radler (1998). In the case of a 
finite fluid body surrounded by free space the electric potential at the boundary has to be 
taken into account. This was already pointed out by Roberts (1967) but not really utilized 
for investigations of specific models. In a recent paper by Dobler and Radler (1998) spatial 
variations of the electrical conductivity were also allowed to occur. However, the case of 
surrounding vacuum had to be excluded in this formulation. 

It is the aim of this paper to re-consider the integral equation formulation of steady 
kinematic dynamo models with finite fluid bodies. This approach may be interesting for 
numerical simulations of dynamos in arbitrary geometry. E.g., realistic simulations of some 
recent laboratory dynamo experiments are still rare due to the fact that for general ge- 
ometry the handling of boundary conditions for the induction equation is not as easy as 
in the spherical case. In calculations concerning the Karlsruhe dynamo experiment it was 
still appropriate to circumvent this problem by virtually embedding the non-spherical (but 
relatively compact) dynamo module into a spherically shaped surrounding of low conduc- 
tivity (Radler et al. 1998). For the Riga dynamo facility with its large aspect ratio another 
method was used (Stefani, Gerbeth and Gailitis 1999). The time-dependent induction equa- 
tion was treated for the dynamo domain with solving, at every time-step, a Laplace equation 
in the outer part with appropriate matching conditions to the inner part. Using this time- 
consuming procedure, a more efficient numerical method which is restricted to the very 
dynamo domain and its boundary seemed highly desirable. In summary, the first reason for 
our interest in the integral equation approach is connected with the search for alternative 
schemes for the efficient and stable numerical treatment of dynamo problems in arbitrary 
geometry. 

The second reason is connected with some new developments concerning inverse prob- 
lems of magnetohydrodynamics (MHD). In two recent papers (Stefani and Gerbeth 1999; 
Stefani and Gerbeth 2000) the problem of reconstructing the velocity field of an electrically 
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conducting fluid from measurements of induced magnetic fields outside the fluid and mea- 
surements of electric potentials at the fluid boundary was addressed. Up to now, this work is 
restricted to small magnetic Reynolds numbers R m , hence an external magnetic field has to 
be applied. The long-term objective is to generalize this inverse problem approach to large 
R m , in particular with regard to laboratory dynamos. Of course, there is a long tradition in 
geophysics concerning inverse problems and (concerning the geodynamo) there is a wide lit- 
erature on reconstructing the tangential velocity at the core-mantle boundary from magnetic 
field observations (see Bloxham 1989 and references therein). However, the frozen-flux ap- 
proximation, which is crucial for that kind of velocity reconstruction, was seriously put into 
question recently (Gubbins and Kelly 1996; Love 1999). For obvious reasons, geophysicists 
never expected that the electric potential at the fluid boundary could also be known from 
measurement. But this situation will be different for laboratory dynamos. Considering the 
huge technical problems for velocity measurements in liquid sodium, an inverse dynamo the- 
ory for velocity determination from measured magnetic fields and electric potentials seems 
to be attractive all the more. For those applications, the integral equation approach with 
its explicit use of the electric potential at the boundary seems to be an appropriate starting 
point. 

The general scheme of the integral equation approach will be represented in section 2. 
The integral equation for the magnetic field is derived from Biot-Savart's law and contains 
a boundary integral over the electric potential. This electric potential, in turn, has to fulfill 
an integral equation over the boundary. It should be pointed out that the incorporation 
of the electric potential at the boundary is well elaborated in the quite different context of 
electro cardiology, electroencephalography, and magnetoencephalography. For dynamos in 
arbitrary geometry, in section 2 also some hints are given concerning the numerical imple- 
mentation of the general scheme. 

In section 3, the general scheme is applied to the case of a spherically symmetric mean- 
field dynamo with an arbitrary radial dependence of a. For that case, we find three coupled 
integral equations for the two defining scalars of the magnetic field and for the electric 
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potential at the boundary. For the special case that a is constant inside the volume we 
re-derive the well-known analytical result obtained by Krause and Steenbeck (1967). 

In section 4 some remarks are made on the generalization of the method to the non- 
steady case and to the case of varying electrical conductivity, and on possible applications 
to inverse dynamo theory. 

II. THE INTEGRAL EQUATION APPROACH 

A. General formulation 

Let us consider a dynamo acting in an electrically conducting non-magnetic fluid which 
occupies a finite domain D with boundary S surrounded by non-conducting space. We 
restrict all the following considerations to the steady case. Then the magnetic field B and 
the electric current density j have to satisfy the Maxwell equations 

VxB = fi j, V-B = (1) 

everywhere with /i being the magnetic permeability of the free space. The electric field E 
has to be irrotational, E = — V</?, with some electric potential if. The current density j is 
assumed to satisfy Ohm's law in the form 

j = a(F - Vy) (2) 

inside D, and it vanishes outside. Here a is the electrical conductivity assumed to be 
constant. F denotes the electromotive force u x B where u is the velocity of the fluid 
motion. In the framework of mean-field electrodynamics (see, e.g., Krause and Radler 1980) 
B, j, tp, and F can be interpreted as mean fields. Then F can be fixed , e.g., to the form 

F = u x B + aB - f3V x B , (3) 

where u denotes now the mean velocity, the term aB describes the a-effect and the term 
(3V x B another effect which can be interpreted by introducing a mean-field conductivity 
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different from a. With a = (3 = we formally return to the case considered before. The 
considerations of this section apply for all F being homogeneous linear functions of B and 
its derivatives independent of their precise form. 

The equations given so far define a problem for B. Together with the requirements that 
there are no surface currents on S and that B vanishes at infinity they allow to determine 
B (apart from a constant factor) if the dependence of F on B is fixed. 

As it is well known equations (1) are equivalent to Biot-Savart's law 

= £ I j(r') x (r - rO dv , 
4:71 J d \r - r'\ 3 

With j according to (2) we obtain 

= ^ r F(r|)x(r r'| ^ w [ y r-s_ ^ 

4:71 Jd |r-r'| 3 4n Js \r - s'\ 3 

with n(s') denoting the outward directed unit vector at the boundary point s' and dS' 
denoting an area element at this point. 

In contrast to the differential equation approach which usually deals with the magnetic 
field only, we have now to deal with both B and the electric potential (p which is, however, 
needed only at the very boundary. 

According to (1) the current density j is source- free everywhere. Therefore we may 
conclude from (2) that 

Aip(r) = V • F(r) (6) 

in D. Since there are no surface currents on S we have there j • n — 0, i.e, the normal 
derivative of the potential on the inner side of S has to satisfy 

dip 



dn 



\ s = F(s) ■ n(s) . (7) 



Using Green's theorem, it can be shown (Courant and Hilbert 1962; Barnard et al. 1967a) 
that 
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47rJ r-r' 47r 7 r - s'\ 

-i- y • j^fpjs dfif (8) 

where p = 1 for points r inside D, p = 1/2 for points r = s on 5 and p = for points r 
outside .D. A solution for ip on 5 can be found by either taking the limit r — > s from outside 
or inside (for the latter case it is important to note that y? is continuous in D + S) or by 
solving the version of (8) for r = s, 

^ W 2tt7 s-r' 2tt7 v ; s - s' 

D S 11 

-i- / ip{s>)n{s>) ■ ^f|p . (9) 

s 

In this context it is of importance that 

hm / y{s')n{s') ■ -^—^ dS' = + I <p(s')n(s') ■ -^-^ dS' (10) 

I J 7» — s \ J \S — S I 

where the upper and lower signs correspond to the approaches to the point s from inside 
or outside D, respectively. The last integrals in (9) and (10) have to be understood in the 
sense of principal values, which are obtained by first integrating over an area S which is S 
diminished by a small tangential disk of radius e around s and taking than the limit e — > 0. 

The two integral equations (5) and (9) provide another complete formulation of the 
problem for B as it was defined above on the basis of differential equations. 



B. Some numerical aspects 

Let us make some formal remarks concerning the numerical implementation of the cou- 
pled system of equations (5) and (9) for dynamo problems in arbitrary geometry. Assume 
certain discretizations of integrals and let us denote all components of B at the grid points 
i by Bi as well as all ip at the boundary grid points m by <p m . The strength of the induction 
effects incorporated in F is scaled by a factor A. Using Einstein's summation convention, 
equation (5) can formally be written as 
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B t = \M tk B k + N im ip m . (11) 

For any chosen discretization, the precise form of the matrices M and N can easily be 
derived from equation (5). We want to point out that M depends on F but N depends 
only on the geometry of the boundary. 

The integral equation (9) for the electric potential can as well be written in matrix 
notation as 

(fi + Ei m (p m = \H ln B n (12) 



or 



Glmtpm = ^H hl B n (13) 



with 



Glm — Sim + Elm ■ (14) 

The matrix G depends only on the geometry of the boundary. Concerning the solution of 
(13) some care is needed as G is singular. This is connected with the fact that (p is only 
determined up to a constant. However, there exist methods to circumvent this problem, 
one of them being the deflation method (see, e.g., Barnard et al. 1967b) where the singular 
matrix G is replaced by an appropriate non-singular matrix G giving a unique inversion. 
Rather than going into the details of the deflation method, let us assume for the moment 
that we have found the inverse of G. Then ip can formally be written as 

if m = \(G- l ) m iH ln B n . (15) 

In a last step this can be inserted into equation (11) to give the matrix eigenvalue equation 

B, t = \(M lk B k + N m (G- l )miH ln )B n . (16) 

Note again that N and G 1 depend only on the geometry of the boundary. In principle, 
when dealing with various u, a and (3 the product matrix iV • G^ 1 must be computed only 



once for a given geometry. 

The preceding considerations provide the framework for numerical computations for finite 
domains of arbitrary shape. 

III. A SPHERICAL MEAN-FIELD DYNAMO MODEL 

We apply our integral equation approach now to a simple spherical mean-field dynamo 
model which is a slightly modified version of a model proposed by Krause and Steenbeck 
(1967) (see also Krause and Radler 1980). To define this model we specify D to be a spherical 
region and put F = aB, i.e., u = and (3 — 0. In contrast to the original model, whose 
advantage is the possibility of an analytical treatment, we consider a no longer as necessarily 
constant but admit it to vary with the radial coordinate r. Starting from equations (5) and 
(8) we will derive three integral equations for three functions of the radial coordinate r, two 
of which define B and the third one tp. In particular, we will re-derive the known analytical 
result for the original model within this new approach. 

We want to point out that our model involves quite a few simplifications, and therefore 
the results have to be considered with care. In particular, an ideal a-effect as given by 
F = aB occurs only with a homogeneous isotropic turbulence and then a has to be con- 
stant. A spatial variation of a requires deviations from this assumption which necessarily 
leads to other contributions to F, which are ignored here. Those deviations from homoge- 
neous isotropic turbulence and additional contributions to F must of course occur near the 
boundary of the fluid body. In the original model of Krause and Steenbeck this neglect led 
to a conflict between the results obtained for the high conductivity limit and a statement 
by Bondi and Gold (1950) according to which in this limit, roughly speaking, any growth 
of the magnetic field in outer space has to be excluded. Actually, this conflict was resolved 
by a consequent treatment of a model taking into account additional contributions to F 
(Radler 1982), and it was shown that nevertheless a dynamo is well possible even in the 
high-conductivity limit but the magnetic field is then completely confined inside the fluid 
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body (Radler and Geppert 1999). 



A. Mathematical preliminaries 

We start with splitting the magnetic field B into poloidal and toroidal parts, B P and 
Bt, and representing them by defining scalars S and T, 

S P = Vxv(^r), B T = v(^r) . (17) 

We refer here to spherical coordinates r, 9, 4> and denote the radius vector by r. The defining 
scalars and the electric potential are expanded in series of spherical harmonics Y tm , 

S(r, M) = E sim(r)Y lm (9, 0), T(r, 0,0) = £ *im(r)lW0, 0), 

p(r,0,0) = £^ m (r)Y Jro (M)- (18) 

l,m 

The Yz m (0, 0) are defined as 

' 2Z + 1 (Z-m)! 



P, m (cos0)e im *, (19) 



\ 47r (Z + m)! 

with P™ being associated Legendre Polynomials. The summation in (18) is over all Z and 
m satisfying I > and \m\ < I. Since, however, terms with I = are without interest in the 
following we restrict all discussions to I > 1. Since S, T and ip are real we have si- m = s^ m 
and analogous relations for ti m and ipi m . The definition (19) implies 

/ rf0 / sin0ci0F; m ,(0,0)y m (0,0) = 8 lv 8 mm , . (20) 

JO JO 

In addition we have 

n Y lm = -1(1 + l)Y lm (21) 

where the operator Q is defined by 

1 d { d j~ \ 1 /" 

+ ^^- (22) 
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From (17) and (18) we obtain with the help of (21) the components of B 

B r (r, M)=E l -^^-s lm (r)Y lm (e, (j>) (23) 

l,m r 

f^\rsmd o<p r dr o9 ) 

B ( r 0A\^y^( tim(r)dY lm (9,(p) 1 ds lm (r)dY lm (e,(f)) \ 
* ' ' r 99 r sin 6 dr d<p J 

Finally we recall the expression for the inverse distance between two points r and r', 

I OO I 1 I 

173^ = 4 -E E ^—^YLio^mUOA) (26) 

where r > denotes the larger of the values r and r', and r< the smaller one. 

In the following we will derive integral equations for the functions Si m (r) and t/ m (r), and 
relations for the coefficients (pi m (R). 

B. The defining scalar of the poloidal part of the magnetic field 

Taking the scalar product of both sides of (5) with the unity vector e r we obtain 



= ^ /■ a(rOB(rO x (r - r>) _ ^ , x . „.,, 

47T Jd |r-r'| 3 An Js \r - s'\ 3 

Hqct f V r > x (a(r')B(r')) r' 



I v - x ^rmrn . e , 1 dr . (27) 

47r |r — r'| r 

In the derivation of the second line of (27) we have expressed e r under the integrals by 
(r — r')/r + (r'/r)e r > and used the fact that n(s') and e r > coincide for r' = s'. Considering 
now 

V r > x (a(r')B(r')) = -B(r') x V r >a{r') + a{r')V r , x B(r') (28) 

in (27) we see that the scalar product of the first term on the right hand side with e r / 
vanishes as the gradient of ct(r') points in r'-direction, too. From (21), (24) and (25) we 
obtain 

(V r , x B(r')) ■ e,, = E W(rO W*', 0') . (29) 



l',m' 
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Taking (27), (28) and (29) together we find 

/ «(0 £ 

I'm' 



^^W'(r') x Y Vm \e>^')—^- dV . (30) 



£ — —s lm {r)Y lm {0,<t>) = --7- 

l,m T r 4f Z'm' 

Expressing the inverse distance according to equation (26) we have to distinguish between 
the cases r > r' and r < r'. After integrating on the right-hand side of (30) over the 
primed angles, multiplying then both sides of (30) with 1^(0,0) and integrating over the 
non-primed angles we obtain the first integral equation of our problem in the form 



Slmir) = Hq(T 



21 + 1 



/ — a {r')t lm {r')dr'+ / — a(r') t lm {r') dr' 
Jo r l J r r 



(31) 



C. The electric potential at the boundary 

For the determination of the potential at the boundary we start from (8) for points r 
outside D. As for the last boundary integral we have 

l - f <p(s')n(s') ■ , r ~ s ' dS' = ±- [ p(s')-?-T-^ dS' 
An J v ' v ' r-s' 3 4tt J rv ' ds' r - s' 



<j (m I'm' ' 

and thus 



hm — / ^s')n(s') ■ ^—^3 dS> = £ — -^(i?)!^, 0) . (33) 

<j r S I lm 

For the evaluation of the volume integral in (8) we use 

V r , • (a(r')B(r')) = (V r >a(r')) • B(r') + a(r')V w • B(r') 

- (34) 

and take B r from (23). In this way we find 
1 r V r , ■ (a(r')B(r')) 1 r da(r>) l'{V + 1) 1 

D D lm 1 1 
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and thus 

1 rV rl -{a{r')B{r'))^ ^ 1(1 + 1) 



lim 

r+s An 



D 



(36) 



Analogously, we obtain for the remaining boundary integral in (8) 

lim - / n(s') ■ a W B W dS> = £ i * Im (i2) 1W(6>, 0) . 



r — s' 



2/ + 1 R 



g || lm 

Evaluating now (8) for r — > s with the help of (33), (36) and (37) we find 

x f R r ' 1 da(r') , , , Z + l 
</W#) = -(/ + 1) y o s im(r ') dr> + — a(i2) s im (i?) . 



(37) 



(38) 



D. The defining scalar of the toroidal part of the magnetic field 



We take now the curl of both sides of (5) thus obtaining 

Vr x B(r) = ^ [v r x (v r x / a ^ B ( r P dV') - V r x / ¥>(*>(«') x ■ ~ ~ ^ ^S" 
47r y |r — r'| y 



r — s 
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• (39) 



Considering first the case r < R we further form on both sides of (39) the scalar product 
with e r . We note that 

e , . ( Vr x V r x / ^g^V) = e„ • ((V,.V r . -A r ) / ^^W') 



<9r 7 



|r — r 

+4?ra (r)B r (r) 



s 



\r — s'\ 



dS' 



(40) 



where we have used the identity A r \r — r'\~ l = —An5(r — r'). The two integrals on the 
second line of (40) were already treated in the last subsection. Concerning the boundary 
integral in (39) over the electric potential we note that 

\ d 2 i 



V r x n(s') X 



r — s 



\r — s 



/|3 



drds' \r — s' 



(41) 



Putting everything together we have 
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(V r x B(r)) ■ e r = fi aa(r)B r (r) + 



/fog; 

47T 



d_ f Mr') B , Tl) 1 ^ d I a(s>)B r (s>) 

dr Jd dr 1 r \r — r'\ dr J \r — s'\ 

ii s i i 



-dS' 



is ' ' ' drds' \r — s' 

Representing the right-hand side according to (29), expressing <p according to (18) and (38), 
using (26) and integrating both sides over the angles we obtain 



dS' 
(42) 



Um{r) = fi a 



l + l r r a da(r') , „ , , / /■« r 1 ^ 1 cta(r') , , , 



21 + 



im»J r Slm{r ] dr ~ a{R) Slm{R) 



From (39) we can also conclude that V x B = for r > R. Without going into the details 
of the proof we note only the consequence 



ti m (r) = for r>R 



(44) 



E. Connection with the differential equation approach 

Notwithstanding the fact that the differential and the integral equation approach are 
equivalent in a general sense it might be instructive to show this equivalence for our spe- 
cial problem. Differentiating equations (31) and (43) two times with respect to the radial 
component we obtain the equations 

d 2 s lm (r) 1(1 + 1) 



dr 2 r 



2 



si m (r) + Hoc a(r) t lm (r) = (45) 



d 2 tim(r) da{r) ds im {r) ^ d 2 s lm {r) 1(1 + 1) 

^ G ~dr dr W<x{r)—^ 2 ^— (ti m {r) - f*r a(r)s lm (r)) = . (46) 

These are (apart from a factor r in the definitions of si m and t[ m ) the same differential 
equations for the considered problem of radially varying a for the steady case as they were 
already derived by Radler (1986). 

The boundary conditions which are used in the differential equation approach can as 
well be derived from (43) and (31). For the case r = R we see that the third term on the 
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right hand side of (43) vanishes identically and that the second and fourth term as well as 
the first and the fifth term are canceling each other. Thus we arrive in a natural way at 



tlm(R) = 



(47) 



which is one of the boundary conditions. The second boundary condition 

ds im(r) ■ I (R) _ n 

1 \ r =R + -Slm{K) — U . 

dr r 



(48) 



can be derived by differentiating equation (31) with respect to the radius and using (47). 



F. Result for the dynamo model of Krause and Steenbeck 



Let us now specify our results to the original model by Krause and Steenbeck, i.e., to 
the case of constant a. From equations (31), (38), and (43) we obtain 



1 



21 + 1 



/ —j-t lm {r')dr' + / —-t lm {r')dr' 
Jo r l J r r 



tim{r) 



jj aa (si m (r) - j^s lm (R)^ for r < R 
for r > R 



(49) 
(50) 
(51) 



The equation (50) is already incorporated in (43) or (51) and is not needed for the de- 
termination of the magnetic field but allows to calculate the electric potential afterwards. 
Introducing the dimensionless variable x = r/R we can rewrite now the two equations (49) 
and (51) for the expansion coefficients of the defining scalars in the form 



sim(x) = /j, aa 



R 2 



21 + 1 



/ r t lm (x') dx' + / —T-tiraix^dx' 

JO X 1 Jx X n 



tlm(x) 



J /iocra (si m (x) - x l+1 si m (R)^j for r < R 

I 



for r > R 



(52) 
(53) 



These equations are solved by 

c ^ (j^^^Ji+i/zil^ovaRx) - 



Slm(x) = 



—x l+1 



21+1 



2!+~i x l Ji-i/2{Ho&®-R) 



Ji-i/i{lJ>QaotR)} for r < R 
for r > R 



(54) 
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tl m {x) 




x l 2 Ji+i/ 2 (ijLocraRx) for r < R 
for r > R 



(55) 



if the condition 



Ji+i/2{^olR) = 



(56) 



is satisfied. Here Q m means an arbitrary constant and the J n are Bessel functions of the first 
kind with the order n. This result can be proved with the help of the known relations 



Apart from a slight difference in the definitions of the defining scalars (by a factor r) the 
result (54)-(56) coincides with that reported by Krause and Radler (1980). 



We have restricted our considerations to the steady case and to a constant electrical 
conductivity of the fluid. It is easy to modify them in such a way that they apply to 
a spatially varying conductivity. Then the electric potential is not only needed at the 
boundary of the fluid but at all places with a non-zero gradient of the conductivity. As 
for the non-steady case, it should be noted that for an infinitely extended fluid Dobler and 
Radler (1998) were able to formulate an integral equation approach under the assumption 
that the time dependence of the magnetic field has the form B(r,t) = exp (jt) B{r). A 
similar procedure should also be possible for time-dependent dynamos in finite domains. 

As already mentioned in the introduction, on of the long-term goals is to develop an in- 
verse dynamo theory. For small R m , using an applied external magnetic field, it was shown 
(Stefani and Gerbeth 2000) that the velocity field can be reconstructed from the measure- 
ments of the induced magnetic field and the induced electric potential at the boundary if 
some regularization of the inverse problem is used. The presented integral formulation is 




(57) 



Jn+l{x) + Jn~l{x) = 2UX 1 J n {x) . 



(58) 



IV. GENERALIZATIONS AND PROSPECTS 
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intended also as a first step to generalize this method to large R m . Still needed is an appro- 
priate formulation of the inverse dynamo theory in terms of an eigenvalue equation combined 
with a least-square adjustment calculus giving an optimal fit of the model to the measured 
quantities. Surely, some kind of regularization (e.g. with some quadratic functionals of the 
velocity or the magnetic field) will be needed for that purpose. Ideas pointing in a similar 
direction are discussed by Love and Gubbins (1996). For those investigations in inverse 
dynamo theory the presented integral formulation of the forward dynamo problem might be 
useful. 
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